Can a changing a explain the Supernovae results? 

John D. Barrow^ and Joao Magueijo^ 
^DAMTP, Centre for Mathematical Sciences, Cambridge University, Wilberforce Rd., Cambridge 

CBS OLA, U.K. 

' ^Theoretical Physics Group, The Blackett Laboratory, Imperial College, Prince Consort Rd., 

o ■ 

(N ■ London, SW7 2BZ, U.K. 

X) 



(N 
(N 

> 

m 
o 

o 



Received ; accepted 



- 2 - 
ABSTRACT 

We show that the Supernovae results, implying evidence for an accelerating 
Universe, may be closely related to the recent discovery of redshift dependence in 
the fine structure constant a. The link is a class of varying speed of light (VSL) 
theories which contain cosmological solutions similar to quintessence. During the 
radiation dominated epoch the cosmological constant A is prevented from dominating 
the Universe by the usual VSL mechanism. In the matter epoch the varying c effects 
switch off, allowing A to eventually surface and lead to an accelerating Universe. By 
the time this happens the residual variations in c imply a changing a at a rate that is 
in agreement with observations. 

Subject headings: Cosmology: theory - observations - early universe 
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1. Introduction 

Two puzzling observations are challenging cosmologists. The Supernovae Cosmology Project 
and the High-z Supernova Search (Perlmutter et al 97, Garnavich et al "9^ , ^chmidt 98 , [Riess 



et al 98|) have extended the reach of the Hubble diagram to high redshift and provided new 



evidence that the expansion of the universe is accelerating. This may imply that there exists 
a significant positive cosmological constant, A. In separate work, the spacings between quasar 
(QSO) absorption lines were examined in Keck I data at medium redshifts, z ~ 1, ( IWebb et al| 



^ and compared with those in the laboratory (see also Drinkwater et al 98 , D amour and Dyson 
96, ^hylakhter 76| , Barrow 87| ). These observations are sensitive to time variations in the value of 
the fine structure constant a = e'^/{hc), (where e is the electron charge, h Planck's constant, and 
c the speed of light), at a rate one million times slower that the expansion rate of the universe. 
Evidence was found for a small variation in the value of a at redshifts ~ 1. This could be produced 
by intrinsic time variation or by some unidentified line-blending effect. In this Letter we assume 
that the variation is intrinsic and show that there may be a link between the observations of 
cosmological acceleration and varying a. 

If A > 0, then cosmology faces a very serious fine tuning problem, and this has motivated 
extensive theoretical work. There is no theoretical motivation for a value of A of currently 
observable magnitude; a value lO^^'^ times smaller that the 'natural' Planck scale of density 
is needed if A becomes important near the present time. Such a small non-zero value of A is 
'unnatural' in the sense that making it zero reduces the symmetry of spacetime. A tentative 
solution is quintessence ( |Zlatev et al 99| ): the idea that A might be a rolling scalar field exhibiting 
very long transients. Here we introduce another explanation. 

There are a variety of possible physical expressions of a changing a. Bekenstein proposed a 



varying e theory ( Bekenstein 82 ). An alternative is the varying speed of light (VSL) theory ( Moffat 



93, Albrecht fc Magueijo 99 , Barrow 99 ) in which varying a is expressed as a variation of the 



speed of light. The choice between these two types of theory transcends experiment, and merely 



reflects theoretical convenience in the choice of units ( Barrow fc Magueijo 981) . The simplest 
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cosmology following from VSL is known to contain an attractor in which A and matter remain at 
fixed density ratios throughout the life of the universe ( parrow fc Magueijo 99| ). Such attractor 
solves the fine tuning problem forced upon us by the supernovae results. Hence there is scope 
for the observed changing a to be related to the observed acceleration of the universe. In this 
Letter we propose a model which leads to good quantitative agreement, given experimental errors, 
between the observations of acceleration and varying a. In Section |2| we examine the construction 
of the Hubble diagram in VSL theories, and the interpretation of varying-a experiments. Then in 
Section |3| we present an example of a VSL model which can jointly explain the supernovae results 
and the Webb et al varying-a results. We conclude with a discussion of some further aspects of 
the model proposed, to be investigated elsewhere. 



2. The VSL Hubble diagram 

The Hubble diagram is a plot of luminosity distance against redshift. The purpose is to map 
the expansion factor a(t), where t is the comoving proper time. Redshifts provide a measurement 
of a at the time of emission. If the objects under observation are "standard candles" (as Type la 
supernovae are assumed to be), their apparent brightness gives their (luminosity) distance, which, 
if we know c, tells us their age. By looking at progressively more distant objects we can therefore 
map the curve a{t). 



We now examine how this construction is affected by a changing c. In Albrecht & Magueijc 
99 we showed that E oc for photons in free flight. We also showed that quantum mechanics 



remains unaffected by a changing c if ?i oc c (in the sense that quantum numbers are adiabatic 
invariants). Then all relativistic energies scale like c^. If for non-relativistic systems h cx 1/c, the 
Rydberg energy = mee^/(2^^) also scales like c^. Hence all absorption lines, ignoring the fine 
structure, scale like c^. When we compare lines from near and far systems we should therefore see 
no effects due to a varying c; the redshift z is still 1 + Zg = ao/ae, where o and e label epochs of 
observation and emission. 



- 5 - 



In order to examine luminosity distances, we need to reassess the concept of standard candles. 
For simplicity let us first treat them as black bodies. Then their temperature scales as T oc 
( Albrecht &: Magueijo J^), their energy density scales as p oc T'^/{hc)^ oc c^, and their emission 
power as P = p/c (X c, implying that standard candles are brighter in the early universe if c < 0. 
However, the power emitted by these candles, in free flight, scales like c; each photon's energy 
scales like c^, its speed like c, and therefore its energy flux like c. The received flux, as a function 
of c, therefore scales like: 



where r is the conformal distance to the emitting object, and the subscripts r and e label received 
and emitted. In an expanding universe we therefore still have 



where PeO is the emitting power of standard candles today. Notice that the above argument is still 
valid if the candles are not black bodies; it depends only on the scaling properties of emitted and 
received power. 

We can now set up the Hubble diagram. Consider the Taylor expansion 



where Hq = do/ao is the Hubble constant, and qo = — doao/do is the decceleration parameter. 
Hence up to second order z = Ho(tQ — t) + (1 + qo/2)HQ{t — to)^, or 




(2) 



a{t) = ao[l + Ho{t - to) - ^qoH^it - t^f + ...] 



(3) 



to - t = —[z - {I + qo/2)z^ + ..]. 



(4) 



-"0 

From (^) we find that the luminosity distance c^l is 




(5) 



The conformal distance to the emitting object is given by r 



= //° c{t)dt/a{t). From (|3|) we have 



that 



r = co[{to -t) + 



1 — n 



Ho{to-tf + ...] 



(6) 



2 
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where we have assumed that locahy c = coa" (that is c = cq[1 + nHQ{t — to) + •••])• Substituting 
(Q) we finally have 

dL = ^[z + \{l-{q, + n))z^ + ..] (7) 

We see that besides the direct effects of VSL upon the expansion rate of the universe, it also 
induces an effective acceleration in the Hubble diagram as an "optical illusion" (we are assuming 
that c decreases in time: n < 0). This is easy to understand. We have seen that VSL introduces 
no intrinsic effects in the redshifting spectral line or in the dimming of standard candles with 
distance and expansion. The only effect VSL induces on the construction of the Hubble diagram 
is that for the same redshift (that is, the same distance into the past) objects are farther away 
from us because light travelled faster in the past. But an excess luminosity distance, for the same 
redshift, is precisely the hallmark of cosmological acceleration. However, we need to consider 
the other experimental input to our work: the [Webb et al 99 results. By measuring the fine 
structure in absorption systems at redshifts z ~ 0(1) we can also map the curve c(t). Since 
c = co[l + nHo^t — to) + ...] we have c = co[l — nz + ...], and so to first order a = ao[l + 2nz + ...]. 
However, the results presented in Webb et al 99| show that n is at most of order 10~^. This 
means that the direct effects of varying c permitted by the QSO absorption system observations 
are far too small to explain the observed acceleration. We need to look at a fully self-consistent 
generalisation of general relativity containing the scope for varying c. 



3. The model 

We start with some general properties of the dynamics of c. Drawing inspiration from dilaton 
theories (like Brans-Dicke gravity) we take ip = log(c/co) as the dynamical field associated with 
c. Indeed, powers of c appear in all coupling constants, which in turn can be written as e'^, where 
(f) is the dilaton. Another theory using a similar dynamical variable is the changing a theory of 

-"^Had we assumed that /i oc c for all systems we would have got instead dL = (co/i^o)[-2 + ^(1 ~ 
(go(l + 4n) + n))z'^], with Hq = (1 — An)HQ. This does not affect any of the conclusions. 
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Bekenstein 82 (which uses log (a)). 



We then endow tp with a dynamics similar to the dilaton. The left-hand side for the tp equation 
should therefore be Oip (in the preferred Lorentz frame - to be identified with the cosmological 
frame). This structure ensures that the propagation equation for -0 is second-order and hyperbolic, 
i.e. propagation is causal. Since VSL breaks Lorentz invariance other expressions would be possible 
- but then the field ^p would propagate non-causally. An example is {g^^ + u^Uu){V fj_ip){Vu'4'), 
where is the tangent vector of the local preferred frame. 

On the other hand one need not choose (as in Brans-Dicke theories) the source term to 
he p — 3p, where p and p are the energy density and pressure of matter respectively. Without 
the requirement of Lorentz invariance other expressions are possible, and using them does not 
conflict with local causality. If T^'^ is the stress-energy tensor we can choose as a source term 
T^^{g^u + u^Uy)] that is changes in c are driven by the matter pressure. We find this choice is the 
one that gives interesting effects. 

For a homogeneous field in an expanding universe we therefore have -0 -|- 3(d/a)?/' = AttGlop/c^ , 
where p is the total pressure of the matter fields and u is a coupling constant (distinct 
from the Brans Dicke coupling constant). The full self-consistent system of equations in a 
matter-plus-radiation universe containing a cosmological constant stress pA = (Ac^)/(87rG) is 
therefore 

+ 3-0 = iTTGu;^, (8) 
a 3 

a ■ 

Py + A-p^ = -2pA0, (9) 
a 

PA = 2/9A0, (10) 

a 

Pm + 3-pm = 0, (11) 



a 

a\2 8-kG 



{pm + Pj + pa), (12) 



,a/ 3 

where subscripts 7 and m denote radiation and matter respectively. We have assumed that the 
sink term ( p!o| ) is reflected in a source term in (^) (and not in (11)). This is due to the fact that this 



term is only significant very early on, when even massive particles behave like radiation. We have 
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ignored curvature terms because in the quasi-lambda dominated solutions we are about to explore 
we know that these are smaller than p\, ( parrow fc Magueijo 99 ). Here, in complete contrast to 
Brans-Dicke theory, the field ifj is only driven by radiation pressure in the dust-dominated era. In 
other words, only conformally invariant forms of matter couple to the field ip. 

In a radiation-dominated universe the behaviour of this system changes at the critical value 
Lo = —4. For u) < —4 we reach a flat pa = attractor as t ^ oo. For —4 < a; < we have 
attractors for which pA and maintain a constant ratio (see [Barrow &: Magueijo 99| ). In Fig. [l| 
we plot a numerical solution to this system, with uj = —4.4 (a 10% tuning below the critical 



value UJ = —4) and n = —2.2 during the radiation epoch. As expected from Barrow &: Magueijo 
99, this forces f^A to drop to zero, while the expansion factor acquires a radiation-dominated 
form, with a oc t^^"^. By the time the matter-dominated epoch is reached, 0,\ is of order 10^^^. 
During the matter epoch, the source term for tp disappears in eq. (|8|), n starts to approach 
zero, Cl\ starts to increase, and the expansion factor takes on the a oc t^/^ dependence of a 
matter-dominated universe. A few expansion times into the matter epoch, becomes of order 1 
and the universe begins accelerating. By the time this happens n is of order 10^^, in agreement 



with the expectations of Webb et al 99 . This type of behaviour can be achieved generically, for 
different initial conditions, with a tuning of to that never needs to be finer than a few percent. 

We can provide an approximate argument explaining why this theory should display 
this type of behaviour and why we need so little fine tuning of uj to explain the supernovae 
experiments. If we neglect changes in c after matter-radiation equality, tgq, we are going to require 
Pkiieq) / p{ieq) ~ ~ lO""*^^. Let c = CQfl'^^*-' , with n = —2 — 5, and n = uj/2 during the radiation 
epoch. We can integrate the conservation equations to give 



n 



PA a^PA n + 2 
with A constant, from which it follows that 



(13) 



PA 



1 + 



S Pi 

2pAi 



25 



(14) 



We see that assymptotically p/pA grows to infinity, if 6 > (the flat pa = attractor of Barrow 
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& Magueijo 99). However the growth is very slow even if 6 is not very small. Our theory 
displays very long transients, and a very slow convergence to its attractor, a property similar to 
quintessence models ( ^latev et al 9£ ). It is therefore possible to achieve pa/p ~ lO^^'^ at the end 



of the radiation epoch, with 6 chosen to be of order 0.1. 



Now, why is the change in c of the right order of magnitude to explain the results of Webb et 
al 99? With a solution of the form c = coa"'-*'* we find that 



n{t) « -—^ (15) 

3(p^ + 2pa) ^ ^ 

valid in the matter dominated era, regardless of the details of the radiation to matter transition. 
With ui ~ —4 we therefore have 

4 2.3 X 10-5 , , 



of the right order of magnitude. The order of magnitude of the index n ~ 10 ^, observed by [Webb 



et al 99, is therefore fixed by the ratio of the radiation and the matter energy densities today. 



4. Discussion 



In this Letter we proposed a theory relating the supernovae results and the observations by 



Webb et al 99. The theory we have proposed is one example within a class whose members exhibit 
similar behaviour. In these theories the gravitational effect of the pressure drives changes in c, 
and these convert the energy density in A into radiation. Thus pA is prevented from dominating 
the universe during the radiation epoch. As the universe cools down, massive particles eventually 
become the source of pressureless matter and create a matter-dominated epoch. In the matter 
epoch the variation in c comes to a halt, with residual effects at z « 1 — 5 at the level observed by 
Webb et al. As the c variation is switched off, the A stress resurfaces, and dominates the universe 
for a few expansion times in the matter-dominated era, in agreement with the supernovae results. 

In a forthcoming publication we shall address other aspects of this theory, beyond the 
scope of this Letter. We mention nucleosynthesis, the location in time of a quantum epoch. 
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and perturbations around the homogeneous solution discussed here (see Barrow & O 'Toole 
99). Nucleosynthesis in particular may provide significant constraints on this class of models. 
However we expect a variation in a to require variations in other couplings if some unification 
exists. Nucleosynthesis involves many competing effects contributed by weak, strong and 
electromagnetic, and gravitational interactions and we do not know how to incorporate all the 
effects self-consistently. Studies of the effects of varying constants coupled by Kaluza-Klein extra 



dimensions have been made by |Kolb et al Sq and [Barrow 87| . The most detailed study to date was 
conducted by Campbell and Olive 95| . 
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Fig. 1. — Evolution oi t x h (where /i = = and JIa in log or linear plots as appropriate. 
The panels on the left (right) describe the radiation (matter) dominated epoch. We have taken 
u = —4.4 (a modest 10% tuning over the critical value u> = —4). In the radiation epoch n = —2.2, 
fl\ slowly drops to zero, and the expansion factor has the usual dependence a oc t^/^. As the 
Universe enters the matter epoch n starts dropping towards zero, a oc t"^/'^, and then starts to 
increase. Eventually n is of order 10~^ and fl\ 0.7. This type of behaviour occurs for a large, 
non finely tuned, region of couplings uj and initial conditions (here Q\ = J7-y = 0.5). 



